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ABSTRACT 

Context. Most multi-planetary systems are characterized by hot-Jupiters close to their central star, moving on eccentric orbits. From 
a dynamical point of view, compact multi-planetary systems form a specific class of the general A'-body problem (where N > 3). 
Moreover, extrasolar planets are found in prograde orbits about their host star, and often in mean motion resonances (MMR). 
Aims. In a first step, we study theoretically a new stabilizing mechanism suitable for compact two-planet systems. This mechanism 
involves counter-revolving orbits forming a retrograde MMR. In a second step, we investigate the feasibility of planetary systems 
hosting counter-revolving planets. Dynamical stability, observations, and formation processes of these systems are analyzed and 
discussed. 

Methods. To characterize the dynamical behavior of multi-dimensional planetary systems, we apply our technique of global dynamics 
analysis based on the MEGNO indicator (Mean Exponential Growth factor of Nearby Orbits) that provides the fine structure of the 
phase space. In a few cases of possible counter-revolving configurations, we carry out new fits to the observations using the Pikaia 
genetic algorithm. A statistical study of the stability in the neighborhood of different observed, planetary-systems is completed using 
a Monte-Carlo method. 

Results. We analyse the observational data for the HD 73526 planetary system and find that counter-revolving configurations may be 
consistent with the observational data. We highlight the fine and characteristic structure of retrograde MMRs. We demonstrate that 
retrograde resonances open a family of stabilizing mechanisms involving new apsidal precession behaviors. 

Conclusions. Considering two possible formation mechanisms (free-floating planet and Slingshot model), we conclude that counter- 
revolving configurations are feasible. 

Key words, celestial mechanics - planetary systems - methods: numerical, observational, statistical 



1. Introduction 

At present, 27 1 extrasolar planets have been detected around 233 
stars (both solar and non-solar typefl Among them, there are 25 
multiple-planet systems : 17 two-planet systems (e.g. HD 82943, 
47 UMa, HD 108874, HD 128311), 6 three-planet systems (e.g. 
V And, HD 69830, Gliese 876, Gliese 581), 1 four-planet system 
(HD 160691) and more recently 1 five-planet system (55 CnC). 
Observations indicate that Mean Motion Resonances (MMR) 
frequently occur for planets of multiple-planet systems : Gliese 
876 (e.g. Rivera et al. 2005), HD 82943 (e.g. Ji et al. 2003, 
Mayor ef al. 2004) and HD 128311 (Vogt et al. 2005) are in 2:1 
MMR, HD 202206 is in 5:1 MMR (Correia et al. 2005), while 
47 UMa is close to a 7:3 (Fischer et al. 2002) or 8:3 commensu- 
rability (Fischer et al. 2003). 

This work is devoted to compact multi-planetary systems, 
characterized by (a) giant Jupiter-like planets found close to their 
central star, and (b) high eccentricities. These two peculiarities 
lead to strong gravitational interactions between the planets and 
may result in an unstable, dynamical behavior. However, we 
observe many such planetary systems suggesting that they are 
stable, and raising the question of why they are stable. From a 
dynamical point of view, compact multi-planetary systems form 
a specific class of the general N-body problem (with N > 3) 
whose analytical solutions are not necessarily known. A stabil- 
ity analysis of planetary systems, using numerical methods to 
explore multi-dimensional parameter space, typically leads to 



stability maps in which rare islands of stability can be identi- 
fied amidst large chaotic zones. The underlying mechanisms for 
these stability zones must be identified. 

In 2002, Kiseleva-Eggleton et al. (2002) showed that the 
cuiTently-published, orbital parameters place the planetary sys- 
tems HD 12661, HD 38529, HD 37124, and HD 160691 
in very different situations from the point of view of dy- 
namical distribution. Since this first study of the com- 
parative stability of multi-planetary systems, many studies 
have been carried out in this direction. The role of the orbital 
mean motion resonances, in particular with a 2:1 ratio, has 
been intensively studied by several research groups (for exam- 
ple Hadjidemetriou 2002; Lee & Peale 2002, 2003; Bois et al. 
2003; Ji et al. 2003; Ferraz-Mello et al. 2005b; Psychoyos & 
Hadjidemetriou 2005; Beauge et al. 2006). As a result, it has 
been discovered that an extrasolar planetary system, even with 
large planetary masses and eccentricities, can be stable if plan- 
etary orbits are close to stable, resonant, and periodic orbits. 
It has also been established (see e.g. Chiang & Murray 2002; 
Lee & Peale 2002; Libert & Henrard 2006) that orbits in a large 
number of compact multi-planet systems, are locked in Apsidal 
Synchronous Precessions (ASP hereafter), i.e. that the apsidal 
lines precess, on average, at the same rateQ. A solution involving 



January, the 11*^^, 2008. |http://exoplanet.eu/catalog.php | 



- For the study of 3-D, full 3-body problems, we introduced the ter- 
minology ASP, for expressing that the apsidal lines precess on average 
in a 3-D space at the same rate: see Bois 2005, Bois et al. 2005. We note 
that in the planar case this phenomenon is also called "apsidal corota- 
tion" (ACR; Beauge et al. 2003). In a number of papers, one may also 



2 



Gayon & Bois: Are retrograde resonances possible in multi-planet systems? 



both MMRs and ASP describes well the stability of eccen- 
tric, compact multi-planetary systems, but may not however be 
unique. We note, for example, that other multi-planetary systems 
have been found to be mainly controlled by secular dynamics (cf. 
Michtchenko et al. 2006; Libert & Henrard 2006; Ji et al. 2007). 
In the present paper, we illustrate theoretically that other mech- 
anisms can in addition provide the stability in multi-planetary 
systems. 

In the case of the HD 73526 system (2: 1 MMR), Tinney 
et al. (2006) found stability over 1 Myr. Based on the analytical 
classification of Hadjidemetriou (2002) established according 
to a hierarchy of masses and eccentricities, this system could 
instead be classified as unstable. Hadjidemetriou's classifica- 
tion may however be too general to disprove the stability found 
by Tinney et al. (2006). Be that as it may, we use the same 
data as Tinney et al. (see Table [1) and our numerical method 
is outlined in the following section. Exploring the stability of 
the HD 73526 system in orbital parameter space, we find large 
chaotic regions. We find that the published data can even be 
described by a chaotic behavior. We note however that Tinney 
et al. (2006) used a diff'erent definition of stability^ Of course, 
we cannot exclude that the observational data were insuffi- 
cient to allow a reliable orbital fit or the fit itself was not 
adequate. On the other hand, it is also possible that the under- 
lying assumption of two prograde orbits is wrong. When placing 
one of the two planets on a retrograde orbit (which forms 
a system with counter-revolving planets), the stability region 
becomes very large. We will show below that this does not imply 
that the orbital fit is consistent with this stability zone. It implies 
that, in the neighborhood of the observational point, we can 
theoretically find stable solutions for counter-revolving config- 
urations. To distinguish between two resonance cases when both 
planets are in prograde orbits, or when one planet is on a retro- 
grade orbit, we call them prograde and retrograde resonances, 
respectively. 

Presently, all known extrasolar planets in multiple systems 
are believed to revolve in the same direction about their corre- 
sponding central star. Most fitted, orbital elements are derived by 
assuming prograde orbits. This is expected according to current 
theories for planetary formation in a circumstellar disk. In order 
to obtain a planet in retrograde resonance, an additional event 
is necessary such as violent, dynamical evolution of the plane- 
tary system, or a capture of the retrograde planet. In our Solar 
System, comets and the planetary satellites of Neptune, Saturn 
and Jupiter are known to have retrograde orbits. It is, therefore, 
important to investigate the stability of exoplanetary systems 
with a retrograde planet in particular if the observations do not 
yield a stable system when assuming all planets on prograde 
orbits. 

In Section|2] we present our method of global dynamics anal- 
ysis. We show that there exists theoretically initial conditions in 
the vicinity of observational data such that stability is only possi- 
ble for a counter-revolving configuration. It raises the question of 
whether such a configuration is consistent with the observational 

find the incorrect expression "apsidal secular resonances" (ASR). ACR 
and ASP are in general not true secular resonances, as highlighted by 
Ferraz-Mello et al. (2005a). 

^ In the paper of Tinney et al. (2006), the claim of stability is ob- 
tained from the dynamical behaviors of the resonant angles related to 
the 2: 1 MMR, rather than by characterizations of quasi-periodicity of 
the orbital solution. Besides, the notion of stability is only presumed to 
be acquired by the simple absence of planet ejection. We use instead the 
usual definition of stability related to quasi-periodicity (see Section |2j 
and suitable for conservative dynamical systems. 



data of a given system (Section|3]l. In SectionlH we focus on the 
statistical occurence of stable solutions related to both prograde 
and retrograde resonances. This statistical approach is applied 
to three systems in 2:1 MMR and two systems in 5:1 MMR. If 
such systems harboring counter-revolving planets exist, we must 
also consider how they form: we discuss this issue in Section |5] 
By analyzing the parameter space in the vicinity of the best-fit 
of the HD 73526 planetary system, we highlight the fine struc- 
ture of the 2:1 retrograde resonance (Sections |6] and |7]), and the 
nature of associated apsidal precessions (SectionlS). In addition, 
we complete an analoguous study for a theoretical system in 5: 1 
retrograde MMR (Section |9]l. 



2. Method 

In order to explore the stability in the parameter space of known 
exoplanetary system in the case of retrograde resonance, we 
use the MEGNO (Mean Exponential Growth factor of Nearby 
Orbits) method proposed by Cincotta & Simo (2000). This 
method provides relevant information on the global dynamics of 
multi-dimensional, Hamiltonian systems and the fine structure 
of their phase space (Cincotta et al. 2002). It simultane- 
ously yields a good estimate of the Lyapunov Characteristic 
Numbers (LCN) with a comparatively small computational 
effort (Cincotta & Giordano 2000). It provides a clear picture 
of resonance structures, location of stable and unstable periodic 
orbits, as well as a measure of hyperbolicity in chaotic domains 
(i.e. the rate of divergence of unstable orbits). Using the 
MEGNO technique, we have built the MIPS (Megno Indicator 
for Planetary Systems) package specifically devoted to study- 
ing multi-dimensional planetary systems and their conditions of 
dynamical stability. We use the property of stability in the 
Poisson sense: stability is related to the preservation of a neigh- 
borhood related to the initial position of the trajectory. Moreover, 
in the Poincare-Lyapunov sense applied to conservative systems, 
when quasi-periodic orbits remain confined within certain limits, 
they are called stable. We note that chaotic, in the Poincare sense 
means that the dynamical behavior is not quasi-periodic (accord- 
ing to the conventional definition used for conservative dynami- 
cal systems) and does not necessarily mean that the system will 
disintegrate during limited period of time. We have already suc- 
cessfully applied the MEGNO technique with the MIPS package 
to the study of dynamical stability of extrasolar planetary 
systems in a series of previous papers (see e.g. Bois et al. 2003, 
2004). In the MIPS package, let us note that the //, c inclination 
parameters refer to the dynamical, orbital-element independent 
of the sin ii line-of-sight inclination factorQ 

By applying the observational data of the HD 73526 plan- 
etary system (see Table [T]0, and scanning the non-determined 
elements, namely the relative inclination (/,. - i^ - it) and 
the Q.r relative longitude of nodes (Q,. - - Q/,), we find two 
main islands of stability, as shown in Fig. [T^. The first (1) is 
obtained for iV e [8°, 97°], and the other (2) for very high rel- 
ative inclinations, namely iV e [173°, 187°]. In this stability 



In the present paper, masses of planets remain untouched whatever 
the mutual inclinations may be. Our reference frame is related to the 
planetary system itself, then dynamically autonomous relative to obser- 
vations. Moreover, scanning the phase space, our stability maps express 
the variations of two explicit parameters, without implicit and external 
relations. 

' with in addition at / = 0, Qi, = fir = 0, it, = 0, and = 1° (because 
of gravitational interactions of the whole 3-body problem, the relative 
inclination iV = ic - ibi^Q is then free to evolve in a 3-D space). 
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Table 1. Orbital parameters of the HD 73526 and HD 160691 planetary systems. Data come from Tinney et al. 2006 and McCarthy 
et al. 2004 respectively. 



Planets 


(Mo) 


mp sm (; 
{Mj) 


P 

(days) 


a 

(AU) 


e 


(1) 
(deg) 


M 

(deg) 


HD 73526 b 


1.08 + 0.05 


2.9 + 0.2 


188.3+0.9 


0.66 + 0.01 


0.19 + 0.05 


203 + 9 


86+13 


HD 73526 c 


2.5 + 0.3 


377.8 + 2.4 


1.05+0.02 


0.14 + 0.09 


13 + 76 


82 + 27 


HD 160691 b 


1.08 + 0.05 


1.67 + 0.11 


645.5 + 3 


1.50 + 0.02 


0.20 + 0.03 


294 + 9 





HD 160691 c 


3.10 + 0.71 


2986 + 30 


4.17 + 0.07 


0.57 + 0.1 


161 + 8 


12.6+11.2 





Fig. 1. Stability maps for the HD 73526 planetary system in the [iV, D.r] parameter space. Panel (a) is plotted with an integration time 
of 5000 yrs while planel (b) is plotted over one order of magnitude longer, namely 50000 years. One observes the disappearance 
of the stable island (1). Only the island (2) survives for longest timescales. Initial conditions come from the best fit of Tinney et al. 
(2006) (see Table[T]i. Relative inclinations and relative longitudes of nodes are defined as follows : iV - ic - it and Qr = Qc - (by 
convention, at f = 0, = 0° and //, = 0°). Resolution of the grid is 1 10x50 for panel (a), 55x50 for panel (b). Black and dark-blue 
colors represent highly stable orbits (<Y> = 2 + 3% and <Y> - 2 + 5% respectively, <Y> being the MEGNO indicator value) 
while warm colors refer to highly unstable orbits (<Y> » 2). 



map, we highlight that stability does not allow coplanar prograde 
orbits. 

The purpose of fast-chaos indicators, and in particular of 
MEGNO, is to predict dynamical behavior over a long timescale 
using short integration times. Our integration times do not 
mean stability times or prediction limits but, using the MEGNO 
indicator, they express the minimal times for knowing trajectory 
future. As a result of the principle of MEGNO, stability is gener- 
ally acquired for timescales far longer than the integration times. 
The ratio of "prediction time" to "integration time" achieved by 
MEGNO, is optimal. 

The MIPS maps presented in this paper were confirmed by 
a second global analysis technique (e.g. Marzari et al. 2006), 
based on Laskar's (1993) Frequency Map Analysis (FMA). The 
FMA method uses the diffusion rates of intrinsic frequencies as a 
measure for stability. The numerical values of these frequencies 
are provided by this method. The lowest intrinsic frequencies 
determine the necessary integration time. It is, therefore, possi- 
ble that the FMA method requires longer integration times than 
MEGNO. 

While the necessary integration time for FMA is provided by 
the lowest intrinsic frequency, one is a priori free to choose the 
time when applying MEGNO. The advantage is that a shorter 
integration time can be used. On the other hand, there is an 
uncertainty about the good choice of the integration time which 
might be too short. We, therefore, produce maps at different 



instant of times and consider the evolution of the most stable 
regions. Fig. [1^ and Fig. [TJi show newly-obtained [iV, Q,] maps 
for the HD 73526 system at 5000 and 50000 years respectively. 
While the first island (1) of Fig. [T^ completely disappears in 
Fig.lTJ), the second one remains always highly stable. Continuing 
in time, island (2) persists. We also find three very small stable 
islands (3), (4) and (5). One of the islands (3) is separated by 
180° in Q.r from the center of the large island (2). The two other 
islands, (4) and (5), are distributed symmetrically with respect 
to the islands (2) and (3). As a consequence, due to the life- 
time and size of each stability zone, the large island (2) with 
ir e [173°, 187°] and Q,. € [173°, 266°] contains the most stable 
orbits (i.e. the least "model dependent" on added perturbations)]^ 
This does not necessarily imply that the two observational 
planets of the HD 73526 system are counter-revolving planets^ 
As a consequence, we study in the following section whether 
such a counter-revolving configuration is consistent with obser- 
vational data. 



Let us note that a relative inclination around 180° is equivalent to a 
planar problem where one planet has a retrograde motion with respect 
to the other. Therefore, considering a scale change of 180° in relative 
inclinations, we will use the notation /J!'"'™ = 1° instead of = 179°, in 
the following. 

' Counter-revolving planets mean that the orbital elements of the two 
planets are orbiting in opposite directions about the central star. 
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3. Observational data fits 

Error-bars in published, orbital elements are significant and a 
fit including new observations may yield quite different orbital 
elements. By using different orbital elements that are consistent 
with observational data, Sandor et al. (2007) found stability for 
the coplanar and prograde case. On the one hand, for lowest 
values of and rms (x^ - 1-57 and rms - 7.9 m.s"' for the 
best dynamical fit of Tinney et al. 2006), we find a weak chaotic 
solution. On the other hand, Sandor et al. (2007) obtain four 
stable solutions but with values ofx^ and rms somewhat higher 
than those of Tinney (x'- e [1.58; 1.87] and rms e [8.04; 8.36]). 
As a consequence, the right astrometric characterization of the 
HD 73526 planetary system still remains open. 

We have performed orbital fits for counter-revolving con- 
figurations using a genetic algorithm (called Pikaia; see 
Charbonneau 1995) based on a fitting methodQ We find stable 
retrograde solutions for values of x^ and rms smaller than the 
prograde fits of Tinney et al. (2006) and Sandor et al. (2007) 
(X^ - 1.257 and rms = 6.34 m.s"'). The radial velocity curve 
of the best stable fit is shown in Fig. |2] It is very similar to the 
radial velocity curve given by Tinney et al. (2006) and Sandor 
et al. (2007). As a consequence, we point out that the possi- 
bility of counter-revolving planets should not be discarded in 
observational-data fits. 

Nevertheless, whatever the directions of motions of the two 
planets are, the x^ values are significantly above the expected 
value of 1.0. More observations would enable better fits to be 
derived. However in these conditions, the possibility may not 
be excluded that the HD 73526 planetary system is a counter- 
revolving system. Anyway, the counter-revolving configuration 
related to the HD 73526 planetary system is consistent with the 
observational data. From a dynamical point of view, counter- 
revolving orbits are all the more plausible because they have 
larger highly stable regions. 




1500 2000 2500 3000 3500 
JD-2441000 



Fig. 2. Dynamical velocity fit with measured velocities of the 
HD 73526 planetary system. The best dynamical fit leading to 
a stable two-planet system is obtained for : nib = 2.4921 My„p, 
ma = 2.5919 Mj„p, Ph = 187.935 days, = 379.795 days, 
ab = 0.6593 AU, a, = 1.0538 AU, = 0.2401, e, = 0.2048, 
= 0°, ic = 180°, Qa = 0°, Q, = 0°, a>b = 184.569°, = 
58.545°, Mb = 97.297°, and M, = 221.361°. The velocity off^set 
is : Vo = -25.201 m.s^\ The corresponding rms residuals are 
6.34 m.s"' while the;t'^ reduced factor is equal to 1.257. 



Orbital fitting process are notably explained in Beauge etal. (2007). 



4. Statistical approach 

By integrating 1000 random systems (according to a Monte- 
Carlo method) within the error-bars proposed by Tinney et al. 
(2006), we obtain statistically more stable solutions for coplanar 
counter-revolving orbits than for prograde ones. For prograde 
coplanar orbits, we find only 17 stable systems while for counter- 
revolving coplanar orbits, we obtain 500 stable systems. 

The occurence of stable counter-revolving systems also 
appears in the neighborhood of other two-planet systems. The 
statistical results for their stability in the prograde case and 
in the counter-revolving one are presented in Table |2] for two 
additional 2: 1 and two 5:1 resonance cases. In all cases, a signif- 
icant number of stable systems in retrograde resonances is found. 
The high statistical occurence of stable retrograde configurations 
justifies the study of such solutions, whether or not they corre- 
spond at present to observational data. On the other hand, we 
are well aware of the cosmogonic problem for obtaining planets 
in retrograde resonances within the frame of current theories of 
planetary formation. During the early dynamical evolution of 
planetary systems, planets may end up on retrograde orbits (as 
discussed in the following section). 



System sources 


Period ratio 


Prograde 
MMR 


Retrograde 
MMR 


HD73526 


2/1 


17 


500 


HD82943 


2/1 


755 


1000 


HD128311 


2/1 


249 


137 


HD 160691 


5/1 





320 


HD202206 


5/1 





631 



Table 2. Statistical results about possibility of stable systems 
to be in retrograde resonance. 1000 random systems have been 
integrated in their errors bars and assuming prograde coplanar 
orbits or retrograde ones. The number of stable systems is indi- 
cated in each case. Data come from Tinney et al. (2006), Mayor 
et al. (2004), Vogt et al. (2005), McCarthy et al. (2004) and 
Correia et al. (2005; table 4) respectively. 



5. Formation of counter-revolving planets 

Up to now, few works have been carried out on the formation 
of highly-inclined or counter-revolving orbits. For instance, 
Thommes & Lissauer (2003) showed that a planetary migration 
leading to resonance capture in 2:1 MMR may cause a signifi- 
cant increase in the mutual inclinations of the planets. However, 
the resulting configuration never seems to exhibit retrograde 
motions. That is why, in this section, we propose two novel 
mechanisms of formation of counter-revolving configurations. 

It has been known for a few years that free-floating 
planetary-mass objects have been located in interstellar space 
(see for instance Zapatero Osorio et al. 2000 and Lucas et al. 
2005). These free-floating planets may interact with planetary 
systems in their host cluster; they are either scattered or captured. 
The recent studies of Varvoglis (2008) show that by integrating 
the trajectories of planet-sized bodies that encounter a coplanar, 
two-body system (a Sun-like star and a Jupiter mass), the proba- 
bility of capture is significant, and almost half of the temporary 
captures are found to be of the counter-revolving type. Although 
captures of free-floating planets remain speculative, this could be 
a feasible mechanism for generating counter-revolving orbits. 

Forming close-in planets by using the slingshot model 
revisited by Nagasawa et al. (2008) is another possibility. 
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Fig. 3. Stability maps in the [ab,eb] and [flced parameter spaces based on the HD 73526 planetary system, taking into account 
initial conditions (C) - 1° and Q,. = 216°). Color scale is the same as in Fig. [1] Resolution of the grid is 50x50. The 

pseudo-observational point is indicated by the intersection of horizontal and vertical lines. V-shape structures correspond to the 2: 1 
retrograde MMR (2:-l MMR). 



Starting from a hierarchical 3-planet system and considering 
a migration mechanism including process of planet-planet scat- 
tering and tidal circularization, the authors show indeed that 
close-in planets may be formed. In a few cases, due to the Kozai 
mechanism (involving exchanges between eccentricities and 
inclinations), one planet enters a retrograde motion. 

Considering (1) these two feasible mechanisms of formation 
of counter-revolving orbits, (2) the dynamical fit obtained in 
Section |3] and (3) the statistical occurence of retrograde solu- 
tions (Section m, we may say that counter-revolving theoretical 
configurations are serious candidates for real systems (that could 
be observed later). In the future, we will study these two forma- 
tion processes more deeply. In the following sections, we focus 
on the specific dynamical behavior of systems harboring 
counter-revolving planets. 

6. Fine structure of the resonance 

Due to the retrograde motion of planet 43, the orbital resonance 
of a given planetary system (e.g. the HD 73526 planetary 
system) is called a 2:1 retrograde MMR (that we also annotate as 
a 2:-l MMR). Studying the fine structure of this MMR provides 
a key to understand the stability of the system. We assume 
initial conditions taken from the stability island (2) of the [iV, Qr] 
map (Fig. [TJ)), that is to say the elements of Table [1] and in 
addition Q,. = 216° and i'';"" = 1°. This set of initial condi- 
tions is noted (C)03- We then obtain the two [at, et] and [gc, ed 
stability maps presented in Fig. [3] The presence of a strong 
MMR generates clear instability zones with a prominent V-shape 
structure in Fig. [3^. We note the narrowness of the V-shapes, 
approximately 0.006 AU wide for the inner orbit (when ei, - 0) 
and 0.0015 AU wide for the outer one (when Cc - 0). As a 
comparison, V-shape structures of the Sun- Jupiter-Saturn system 
are five times and twice as large respectively (for Jupiter and 
Saturn). We note in addition how the pseudo-observational point 
lies at the edge of the V-shapes (Fig.[3]l. 



' By convention, the orbital plane of the planet b is chosen as a refer- 
ence plane. The planet b is supposed moving in the prograde direction 
while the planet c in the retrograde direction. Results remain analogu- 
ous with the reverse assumption (b retrograde and c prograde). 

For n, = 216° and ff"" = 1°, the value of obtained with the 
velocity offset Vo = -38 m.i"' is 2.44. The rms residuals are 12.31. 



7. Evidence for a retrograde resonance 

Maps for extrasolar-planet systems, with stability regulated by a 
prograde MMR, are characterized by small islands of stability 
(or linear stable strips in [flfc,fl(] maps) inside large zones of 
instability. By contrast, in the case of retrograde resonance and 
when assuming the initial conditions (C), we detect a dense, 
stable regime in a series of maps (e.g. [ah, flf]), except for one 
unstable zone related to the MMR. In several cases of resonant 
prograde systems, it has been shown that planets on highly stable 
orbits may avoid close approaches due to their adequate posi- 
tions over their orbits and apsidal line locking (see Bois et al. 
2003). This mechanism of stability is not lost during the dynam- 
ical evolution of the system when the apsidal lines on average 
precess at the same rate (i.e. the ASP phenomenon). Without 
such a protection mechanism of 2:1 MMR combined to an ASP, 
disturbing close approaches between the planets are theoretically 
possible in various planet positions. In the case of the 2:1 retro- 
grade MMR, one planet being retrograde, orbital motions occur 
then in opposite directions. As a consequence, the length of time 
that planets spend in conjunction, is much shorter for counter- 
revolving orbits than for prograde ones. This could explain the 
narrowness of the V-shapes in Fig.|3] 

In Gayon & Bois (2008), we showed that, in cases of very 
compact planetary systems obtained by a scale reduction of a 
given observed system, "retrograde" stable islands survive, in 
contrast to "prograde" ones that disappear This scale reduction 
and the V-shapes of Fig. [3]illustrate the efficiency of retrograde 
MMRs for providing stability. In addition, this MMR mecha- 
nism is coupled to specific behaviors of the apsidal lines, as 
shown in the following section. 

8. A new meclianism of apsidal precession at 
retrograde resonance 

In the case of the 2:1 retrograde MMR (2:-l MMR), the expres- 
sions for the resonance angles 6i and 62, and the ASP angle ^3 
are : 

61 - -Ah - 2Ac + 3(jJh 
02 - -Ah - 2Ac + 3We 

03-3 {(bh - cbc) 
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where A is the mean longitude, and co is the apsidal longitude, 
defined hy A - M + 01 and oj - Q. + (l) respectively, for prograde 
orbital motion of a planet; for retrograde motion, these variables 
are defined to be: A = -M+oj and w = Q-w. The general expres- 
sion for the relative apsidal longitude is: Aoj = w/, - ojc = G^/q, 
where q is the order of the resonance. A more thorough investi- 
gation of retrograde resonances, using an analytical approach, is 
in preparation (Gayon et al. 2008). 
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Fig. 4. Stability map in the [a)h,tOc] parameter space based on 
the HD 73526 planetary system taking into account initial condi- 
tions (C) (/;'■"■'' = 1° and = 216°). The pseudo-observational 
point is inside the stability linear strip. Color scale is the same 
as in Fig.[T] Resolution of the grid is 50x50. 



In Fig. |4]we plot the [a)i,,CL)c] parameter space that shows 
a stable linear strip, in dark-blue, including the "pseudo- 
observational" point. We learn that stable solutions are possible 
only when aii, and ate precess, on average, at the same rate. 
The stabilizing mechanism of the system involves a synchronous 
precession of the apsidal lines. The two longitudes of periastron 
do not precess however in the same direction. The outer orbit 
is aff'ected by a retrograde precession {-cbc) relative to the inner 
orbit's precession (a>h) (see Fig.|5}; vilSh). Writing the longitudes 
of periastron as directed angles, we find that the relative apsidal 
longitude Ad> neither circulates nor librates clearly, as shown in 
Fig.|5h. Aw presents a strange motion composed of : (1) a phase 
of prograde circulation including librations with amplitudes of 
+8°, then following a sharp reversal of circulation direction, 
(2) a phase of fast retrograde circulation, until a second, sharp 
reversal. These two phases alternate successively according to 
alternations (or rocking) of 180°, which correspond to a sort 
of cuspH We note that tbc, in contrast to the case for ajt, does 
not uniformly circulate but presents retrograde circulation phases 
interrupted with short libration intervals ( Fig. ^ vs [Sj^)- In 
spite of the opposite directions of their precession, both orbits 
precess, on average, at the same rate. As a consequence, the plan- 
etary system is affected by an apsidal synchronous precession. 
Considering the alternating behavior of the Aoj angle or the 
unusual presence of cusp in the Ad> behavior, we refers to this 
new stabilizing factor as either an alternating ASP, or a rocking 
ASP (or RASP). We have produced movies illustrating the 



" At phase transitions, a scattering of dots appears. When the outer- 
orbit eccentricity goes to zero, the ojc angle that depends on the 
(a - r)/ae ratio is not defined. 
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Fig. 5. Time variation of the Aw angle (a), cjh (b), and oJc (c) for 
initial conditions (C). In panel (b), the slope of ojh is positive 
while the one of ojc (panel (c)) is negative. Panel (a) expresses 
the behavior of the Aoj combination of ojb and oJc (see text). 



mechanisms introduced in the present paper, in particular the 
phenomenon of RASPp^ 

Relations between the eccentricity of the inner orbit and the 
Aw angle exist. In [Aw, et,] parameter space, we can see from 
Fig- how the 2:1 retrograde MMR spreads out its resources: 

(1) Inside the 2:-l MMR (i.e. inside the [a, e] V-shape), both 
apsidal lines on average precess at the same rate and in the same 
prograde direction. The 2:1 retrograde MMR is then combined 
with a uniformly prograc/e ASP (island (1)). 

(2) Close to the 2:-l MMR (i.e. outside but close to the 
[a, e] V-shape), both apsidal lines on average precess at the same 
rate but in opposite directions. The 2: 1 retrograde near-MMR is 
combined with the mechanism of rocking ASP (island (2)). 

(3) The division between these two islands is related to 
the degree of closeness to the 2:1 retrograde MMR. We high- 
light the fact that, for a long timescale, the third island in the 
[Aw, e;,] map, is proved to be a chaotic zone (island (3) where 
Aw e [80°, 280°]). 



Movies are downloadable from: 
|http://www.oca.eu/gayon/Extrasolar/Retro_MMR/movies.htnil| 
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Fig. 6. Stability map in the [Aw, e/,] parameter space based on 
the HD 73526 planetary system, taking into account initial con- 
ditions (C) (Table [U with in addition = 1°). The pseudo- 
observational point is inside the stability large island (2) char- 
acterized by a 2:1 retrograde near-resonance and an alternating 
ASP. Color scale is the same as in Fig.[T] Resolution of the grid 
is 50x50. 

We note that a mechanism of stability involving an ASP may 
persist far from the MMR in the prograde case, while it disap- 
pears for short distance to the MMR in the counter-revolving 
case. Hence, moving away from the [a,e] V-shape of the 2:-l 
MMR (Fig.lSj, we find that both apsidal lines precess in opposite 
directions but at different rates. 

By studying the parameter-space in the vicinity of the best- 
fit of the HD 73526 planetary-system and searching for stable 
configurations with similar values of nns, a new, theoretical 
mechanism of stability has been discovered. It is characterized 
by a 2:1 retrograde, near-MMR combined to a rocking ASP. 
Such a stability also allows a large range of eccentricities (see 
for instance planet b in Fig.|6]l. Such a mechanism is particularly 
robust on large timescales. It is why such mechanisms involving 
such resources of the 2:-l orbital resonance could prove to be 
relatively generic and suitable for the stability of a class of 
compact multi-planetary systems where other solutions of the 
3-body problem are not possible. 
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Fig. 7. Stability map in the [Aw, eh] parameter space based on the 
HD 160691 planetary system (data from McCarthy et al. (2004) 
shifted to the 5: 1 reti'ograde MMR (i.e. at - 1 .44 AU), see Table 
[T]with in addition i'^.'^"'" - 1°). Islands (1) and (3) are character- 
ized by a 5:-l MMR and an ASP with apsidal alignement while 
islands (2) and (4) by an 5:- 1 MMR and a circulation of the Aoj 
variable. Color scale is the same as in Fig. [1] Resolution of the 
grid is 50x50. 



in 5: 1 retrograde MMR, as for the HD 73526 planetary system in 
2: 1 retrograde MMR, we observe the distribution of the different 
apsidal behaviors (see Fig.|7]i. Islands (1) and (3) are character- 
ized by an ASP with an apsidal alignement. More precisely, both 
longitudes of periastron (w^ and lJc) on average precess at the 
same rate, both in the retrograde direction, while the Aw variable 
librates about 0°: it is a uniformly retrograde ASP. Within 
islands (2) and (4), ibb and w^ precess in opposite directions but 
according to different rates: Aw circulates. Hence, no apsidal line 
locking is required for stability. We find fine, V-shape structures 
in [a,e] maps corresponding to the 5:1 retrograde MMR. The 
width is 0.02 AU for the inner planet (when eb - 0) and only 
0.002 AU for the outer one (when ec - 0). Consequently, the 
dynamical study of this case derived from the HD 160691 plan- 
etary system allows us to find another theoretical possibility of 
stability involving resources of a retrograde MMR. 



9. The 5:1 retrograde MMR 

The mechanism of stability involving both a retrograde MMR 
and an ASP is also found for the 5:1 orbital period ratio by scan- 
ning initial conditions in the vicinity of the HD 160691 planetary 
system (McCarthy et al. 2004, planets b and c, see Table[TJ. By 
detecting a fourth planet, we note that new observations of this 
system have modified the orbital structure of the entire system 
(Pepe et al. 2007). The observations of Pepe et al. (2007) show 
furthermore that a new coplanar fit of prograde orbits and a new 
fourth planet seems to solve the problem without the need for 
retrograde resonance. We note that in the vicinity of the best 
fit solution of McCarthy et al. (2004), it is possible to find an 
example of 5:1 retrograde resonance. We consider this possibil- 
ity as an academic investigation of this order of MMR. 

The important point is that this 3-body system is completely 
unstable for prograde orbits. Nevertheless, by scanning the non- 
determined parameter space, only one island of stability ex- 
ists (for /, e [156°, 204°] and e [7°, 180°]) and after an 
analysis we highlight the retrograde motion of the outer planet. 
Moreover, scanning the [Aw, e/,] parameter space of the system 



10. Conclusion 

We have found novel mechanisms giving rise to stability that 
could be suitable for a class of compact planetary systems. Such 
mechanisms involve counter-revolving orbits forming a retro- 
grade MMR occuring in a quasi-identical plane. High statistical 
occurence of stable counter-revolving orbits is found. Our study 
of retrograde MMRs indicates the large stability domains and 
the specific behaviors of the precession and resonant angles. 
We propose that these large stability domains are caused by 
close approaches much faster and shorter for counter-revolving 
configurations than for the prograde ones. Scanning the 
HD 73526 planetary system, we find evidence for a new type 
of apsidal precession (the rocking ASP). We find that the 
difference between the longitudes of periastron reveals a specific 
alternation mode at retrograde resonances. We emphasize that 
the counter-revolving configuration studied for the HD 73526 
planetary system is consistent with the observational data. Free- 
floating planets or the Slingshot model might explain the origin 
of such counter-revolving systems. 
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